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Summary. — We demonstrate that the apparently ad hoc parametrization of the 
particle production spectra discussed in the literature and used in description of 
cosmic ray data can be derived from the information theory approach to multiparticle 
production processes. In particular, the violation of the Feynman scaling law can 
be interpreted as a manifestation of nonextensivity of the production processes. 

PACS 96.40 - Cosmic rays. 

PACS 01.30.Cc - Conference proceedings. 
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1. Introduction 

The shape of the x = E/Eq spectra of secondaries is of great importance in all 
investigations concerning developments of cosmic ray cascades [1] (cf. also [2]). The 
crucial problem of practical importance is the existence or nonexistence of the Feynman 
scalling, which says that x-spectra of secondaries are energy independent. Because in 
cosmic ray applications one is sensitive essentially to large x region (of, say, x > 0.1), 
the commonly used formula [1] 



(1) 
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stresses this fact by power-like form with exponent obtained from the x ^ 1 limit of 
fragmentation data [1] . Here a' is parameter responsible for the Feynman scaling violation 
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whereas D and a are additional parameters obtained from the fit to data - aU are energy 
dependent (cf. [f] for details). 

We shall demonstrate that (1) emerges in natural way from the information theory ap- 
proach to hadronization processes presented in [3] and extended to allow for the possible 
nonextensivity of the hadronization process by applying Tsallis entropy [4]. (See [5] for 
details of notion of nonextensivity and [6] for its application to high energy reactions.). 

2. Particle production spectra from information theory 

In order to describe the hadronization process in information theory one uses [3] 
the least biased and most plausible single particle distribution f{y) = jj^ (where y 
is rapidity) resulting from hadronization process in which a mass M hadronizes into N 
secondaries of mean transverse mass fir = vA?+(PTp^ each (for simplicity we con- 
sider only one-dimensional hadronization with limited transverse momenta which can, 
however, depend on M and N). This is done maximalizing the Shannon (or Boltzmann- 
Gibbs) information entropy, S = — J dy f{y) ln/(y), under constraints of normalization 
( J dyf{y) = 1 ) and energy conservation {J dyfiT cosh yf{y) = M/N ), which leads to 
the following (extensive) distribution function cf. [3] for details): 

(2) /(y) = ^(^^ cxp [-/3(M, N) ■ iiT coshy] . 

Here Z{M, N) comes from the normalization of f{y) whereas the Lagrange multiplier 
/3(M, TV) is to be calculated from the energy conservation constraint. Notice that there 
is no free parameter here. This should be contrasted with the popular use of eq. (2) as a 
" thermodynamical parametrizations" with inverse "temperature" 1/ (3 — T beeing a free, 
positively defined parameter. In [3] it was shown that, in a wide range of energies and 
multiplicities, (3{M,N) ~ ^ = -i=f . 

However, the comparison with data cannot be done in such model independent way 
because of the fluctuations of inelasticity K (given by inelasticity distribution x(-^) [7]) 
causing fluctuations of M and because of fluctuations of number of particles N — N{M) 
produced from a mass M given by multiparticle distribution P{N] M) [8] (its actual 
shape accounts here for the possible multisource structure of the production process). 
One has therefore 



(3) 



^= / dKx{K)Y,P{N;K^s)N-—^^^—e^^[-P{K^s-N)iiTCOshy\ 
ay Jo ^ Z{Ky/s,N) 



where (3{Ky/s, N) is still calculated from the energy conservation constraint but now for 
a given event, i.e., it is a fluctuating quantity with fluctuations resembling gamma-like 
distribution (and characterised by mean value /3 and normalised variance LOfi). 

It has been shown in [9, 6] that such fluctuations of the parameter of the exponential 
distribution convert it to a power-like distribution: exp(— x/xq) — > exp (— x/xp) = (1 — 
x/xo/a)" with a = il/uip in our case where x — -^coshy and xq is a parameter 

connected with the mean value of /3 ^ -^ . 

This is precisely the form corresponding to distribution obtained using non-extensive 
Tsallis (/-entropy [4, 5], instead of Shannon one (and reproducing it in the limit q —^ 1), 
S, - -(l-/dy[/(2/)]«)/(l-g) ^ 5g=i = - J dyfiy)lnfiy), together with the modified 
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constraint equation J dyiJ^T cosh y[f{y)]'^ ~ M/N. Using it one gets instead of (2) its 
nonextensive version: 



(4) 



fM = 



Zq{M,N) 



[1 - {l-q)|3q{M,N)■^lTC0s\iy] — 



Nonextensivity means that the entropy of the composition [A + B) of two independent 



systems A and B is equal to [5] S\ 



[A+B) 



= 5, 



{A) 



s, 



(B) 



{l — q)Sq -Sq proceeding to 



its usual additive form only in the q — 1 limit. It arises whenever in the system one en- 
conters long-range correlations, memory effects or fractal structure of the corresponding 
space-time or phase space. Such situation is expected to occur also in the hadronization 
processes [6]. 



3. Comparison with experimental data 




Fig. 1. Comparison of eq. (5) with UA5 [10] data for g = 1 and q = 072. 



We make therefore the following conjecture: the convolution (3) can, in practice, be 
replaced by the following simple one parameter formula of the type of eq. (1) which 
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should be used to describe the existing data of [10, 11, 12]: 



(5) 



dy Zq 



1 - {1 - q)f3giyG, -{N{s)))^lT coshy 



The only free parameter (characterising the strenght of fluctuations in the system ac- 
cording to the previous discussion) is the nonextensivity parameter q. Here Zq = 
J dy expq {—PqiiT cosh y) and {N{s)) is the mean (single non-diffractive) charged mul- 
tiplicity at given energy y/s. Notice that because (3q is calculated from the energy con- 
servation constraint which involves all produced particles (i.e., total multiplicity) it is 
calculated here for |(A^(s)) = {Ntotai{s)) particles. For the same reason care must be 
taken when one addresses data at y/s = 630 GeV as part of them is for charged and 
part for neutral particles only. In both cases the /3q must be the same (calculated for 
total multiplicity at given energy) whereas multiplicity in front of the formula has to be 
chosen accordingly to the actual situation. 

In Fig. 1 we show our results both for g = 1 and our best fit with q = 0.72. In all calcu- 
lations the experimentally observed variation of fir with energy has also been accounted 
for by using the following simple interpolation formula: fir = 0.3 + 0.044 ln(\/s/20.) 
GeV. The results are reasonable, especially for 53 and 200 GeV. For higher energies our 
distributions start to be broader than data and this cannot be improved by changing q 
as diminishing its value in order to make distributions narrower will spoil the agreement 
with data for small rapidities. It turns out that UA7 and P238 data cannot be fitted 
together with UA5 data, as they demand a slightly bigger value oi q = 0.85, cf. Fig. 2. 
Notice that P238 data are for charged and UA7 data for neutral particles therefore they 
must be described with, respectively, ^{Ntotai) and ^{Ntotai) in eq. (5), this leads to 
differences clearly seen in Fig. 2. 



4.0 
3.5 
3.0 
2.5 

:§• 2.0 
z 

■^ 1.5 
1.0 
0.5 
0.0 




ir*i^iiir4ir^ 



• P238 + UA7 @ 630 GeV 

(P238) q = 1.0 

(P238) q = 0.72 

(P238) q = 0.85 
-^;^(UA7) q = 0.85 






















q = 0.7 


0.18- 




a 


q=1.0 
q=1.3 


0^„aaa«aaaa»a,,„^ 


- 


0.15- 








0.12- 




J 


ia*"^ ^^ 


\ 


0.09- 




a'* 




a 




i 







% 


0.06- 


» 







a 




a 







% 

O a- 


0.03- 


a 







0.00- 






V 



|y| 



Fig. 2. Comparison of eq. (5) with P238 
[11] and UA7 [12] data for rapidity distri- 
butions. 



Fig. 3. The examples of distribution (4) 
for mass M = 100 GeV hadronizing into 
N = 20 secondaries of (transverse) mass 
lj,rp = 0.4 GeV each. 
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4. Summary and conclusions 

Our approach consists in noticing the striking similarity of eq. (4) to the eq. (1) of [1] 
{dy = dx/x). Fig. 3 shows the characteristic features of fq{y) for the case of hadroniza- 
tion of firebaU of mass M = 100 GeV into A^ = 20 particles of transverse mass /iy — 0.4 
GeV. Notice that for q < 1 one obtains an increase of particle densities in the central 
region connected with its decrease at the edges of rapidity range - clearly resembling 
the characteristic pattern of Feynman scaling violation. We have replaced therefore the 
"exact" formula, (3) by a one parameter fit represented by (5) where q summarizes the 
action of the averaging over the fluctuations caused by the initial conditions of reaction 
represented by the inelasticity and multiplicity distributions. We flnd it very amazing 
that such simple approach coincides practically with empirical formula (1) and with only 
one parameter q describes fairly well all data. Notice that q — 0.72 is not very far from 
q = 0.75, which gives power 1/(1 — q) = 4 in (1). Notice also that our general formula 
allows for description of small x region as well. It would certainly be interesting to con- 
nect q directly to the parameters describing inelasticity and multiplicity distributions or 
to descriptions of leading particles, for example of the type of that presented in [13]. 
Our results seem to indicate that most probably parameter q should be x-dependent re- 
flecting different character of fluctuations of the quantity N/M [9] in the central (mostly 
pionization) and fragmentation regions. This problem will be addressed elsewhere. 
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